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ABSTRACT 



Context. This is the third paper in a series devoted to the Wide-field Nearby Galaxy-cluster Survey (WINGS). WINGS is a long-term 
project aimed at gathering wide-field, multiband imaging and spectroscopy of galaxies in a complete sample of 77 X-ray selected, 
nearby clusters (0.04 <z< 0.07) located far from the galactic plane (\b\ > 20°). The main goal of this project is to establish a local 
reference sample for evolutionary studies of galaxies and galaxy clusters. 

Aims. This paper presents the near-infrared (J,K) photometric catalogs of 28 clusters of the WINGS sample and describes the proce- 
dures followed to construct them. 

Methods. The raw data has been reduced at CASU and special care has been devoted to the final coadding, drizzling technique, 
astrometric solution, and magnitude calibration for the WFCAM pipeline-processed data. We constructed the photometric catalogs 
based on the final calibrated, coadded mosaics (« 0.79 deg") in J (19 clusters) and K (27 clusters) bands. A customized interactive 
pipeline was used to clean the catalogs and to make mock images for photometric errors and completeness estimates. 
Results. We provide deep near-infrared photometric catalogs** (90% complete in detection rate at total magnitudes J « 20.5, K « 19.4, 
and in classification rate at J ~ 19.5 and K as 18.5), giving positions, geometrical parameters, total and aperture magnitudes for all 
detected sources. For each field we classify the detected sources as stars, galaxies, and objects of "unknown" nature. 

Key words. Surveys - Galaxies : Clusters : General - Catalogs 



1. Introduction 

Clusters of galaxies are important for studies of galaxy forma- 
tion and evolution, because they contain a volume-limited popu- 
lation of galaxies observed at the same cosmic epoch. Galaxies 
in clusters are known to follow tight color-magnitude relations, 
which appear to be universal and to have very small intrinsic 
scat ter to the highest redshifts yet observed (see, amongst oth- 
ers. iBlakeslee et al.ll200l iHolden et all 120041 iMei et all 12005 



|Pe Lucia et al. 2007 ). Together with the conventional interpre- 
tation of the c olor-magnitude rela tion as a mass-metallicity cor- 
relation (i.e.. iTrager et al.l 120001) . this implies that the major- 
ity of the stellar populations in early-type cluster galaxies were 
formed via rapid dissipative starbursts at z > 2. Fundamental- 
plane studies of high-redshift cluster galaxies also support 
this conclusion, at least for the more massive objects (i.e., 



Send offprint requests to: Tiziano Valentinuzzi, e-mail: 
tiziano . valentinuzzi@unipd. it 

* Based on observations taken at the United Kingdom Infra-Red 
Telescope, operated by the Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the U.K. 



Ivan Dokkum & Stanfordl2003tlHolden et alj2005l) . although the 
low mass galaxies seem to have undergone m ore extended star 
formation histories (z.e.. |Poggianti et aTkOQll) . 

Theoretically, the existence of such massive and old galaxies 
at high redshift should represent a strong challenge to models 
where galaxies are assembled hierarchically, from a sequence of 
major mergers at progressively lower redshifts. It is not possible, 
however, to exclude by spectrophotometry alone, these galax- 
ies being assembled from sub-units whose star formation has al- 
ready ceased, but which are not accreted until later (the so-called 
'dry' mergers). This is assumed to be the main channel by which 
spheroids grow at z< 1 in the hierarchical scenario. 

On the other hand, if galaxies are formed via mergers, we 
should observe a steady decrease in the mean stellar mass in 
galaxies as we go to earlier cosmic times and the most mas- 
sive m embers of the merger tree would branch into even smaller 
twigs dDe Lucia et alj |2006[) . While it is generally difficult to 
measure galaxy masses, the K-band lum inosity function is be- 
lieved to provide an adequate surrogate dKauffmann & Charlotl 
1 19981) . Indeed, it is known that the rest frame H or K luminos- 
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ity of galaxies is seen to correlate well with stellar and even 
dynamical mass for loc al and high-redshift galaxies (see, i.e., 
iKodama & Bowerll2003l) . 

To reach higher and higher redshifts it is surely necessary 
to adequately constrain galaxy formation models; on the other 
hand, precise knowledge of the properties of clusters and cluster 
galaxies in the local Universe is necessary, as a benchmark for 
higher redshift studies. Th e Wide-field Nearby Galaxy-clusters 
Survey (jFasano et al .120061 hereafter Paper-I) is an answer to this 
need. WINGSQ is a long-term multiwavelength survey specially 
designed to provide the first robust characterization of both the 
photometric and spectroscopic properties of galaxies in nearby 
clusters, as well as determine the variations in these properties. 

The survey core, based on opti cal B,V imaging of 77 nearby 
((z) » 0.05) galaxy-clusters (see, IVarela et all [2009, hereafter 
Paper-II), has been complemented by several ancillary projects: 

(i) a spectroscopic follow up of a subsample of 5 1 clusters, ob- 
tained with the spectrographs WYFFOS @ WHT and 2dF@ AAT; 

(ii) near-infrared (J,K) imaging of a subsample of 28 clusters 
obtained with WFC AM @ UKIRT, presented here; (iii) U broad- 
and H a narrow -band imaging of subsamples of WINGS clusters, 
obtained with wide-field cameras at different telescopes (INT, 
LBT, Bok). The observations and data redu ction for the first 
two follow-up p rojects have been completed (ICava et al.ll2.Q09l : 
iFritz et al1l2009l in preparation; this paper), while the observa- 
tions for the H a and U-band surveys are still ongoing. 

The near infrared section of the WINGS survey ( WINGS - 
NIR) is a collection of ^0.79 deg 2 mosaics in (J,K) bands of 28 
nearby clusters (17 clusters have been observed in both bands). 

WINGS-NIR is by far the largest survey of nearby galaxy 
clusters as far as the area coverage is concerned. In fact, in 
this redshift range, only individual clusters or small cluste r sam- 
ples have been studied in the literature up to no w (e.g., iPahrd 
U999tlGavazzi et al.lll990tlDe Propris et al.ll2003l) . The WINGS 
survey of near-infrared data consists of nearly one million de- 
tected sources, with 150,000 and 500,000 well classified stars 
and galaxies, respectively. 

In section II of this paper we describe observations and data 
reduction techniques, including a brief presentation of both the 
instrumentation and the software pipeline used to create the final 
coadded mosaics. Section III presents the step-by-step procedure 
of the catalog production pipeline, including an extensive discus- 
sion about star/galaxy classification and interactive cleaning of 
the catalogs. Section IV deals with data quality assessment and 
overall properties of the catalogs, like completeness, astrometric 
and photometric accuracy and precision. In the last section we 
present a brief summary of the main features of WINGS-NIR 
survey and catalogs. 

Throughout this paper we will use a cosmological model 
with parameters: Hq = 70 kms -1 Mpc~', fl^ = 0.3 and £2a 
0.7 



2. Observations and Data Reduction 

The near-infrared data have been collected at the UKIRT tele- 
scope using the WFCAM instrument (see, ICasali et"aTI 120071) 
during four observing semesters, from April 2005 to April 2007 
(see Tab[TJ. The images where obtained in the J and K broad- 
bands of the Mauna Kea photometric systerrQ dTokunaga et al.l 
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Fig. 1. Distribution of WINGS-NIR clusters (red full line) com- 
pared with the distribution of all the WINGS survey clusters 
(black dotted line). The dotted diagram is slightly shifted for 
ease of view. 



120021) . The original plan of the WINGS-NIR survey was to im- 
age at least all the WINGS clusters with spectroscopy and visible 
imaging by UKIRT: unfortunately bad weather (see Tab[T|) lim- 
ited the sample to only 28 clusters. However, as it can be seen in 
Figfflthe NIR subsample has good coverage of X-ray luminosi- 
ties and redshifts when compared with the total WINGS sample, 
while it is slightly biased towards the low cluster velocity dis- 
persion regime. 

WFCAM is an assembly of four Rockwell Hawaii-II 
2048x2048 18jt/m-pixel array detectors with 0'.' 4 pixel size and 
a total field of view of 0.21 deg 2 . The four detectors are spaced 
at 94% of their width, at the corners of a field of about 0.4 deg 2 . 
More details about the ca mera, imaging pr operties, defects and 
drawbacks can be found in lDve et"alT(l2006l) . Figf3]shows the an- 
gular vs. physical median FWHMs for J- (blue open circles; dashed 
line) and K-band (red full circles; full line). It is seen that more 
than 50% of our images have seeing below l'.'O, and only three 
out of 47 co-added mosaics do not meet the quality requirements 
of the WINGS survey regarding physical size resolution (dot- 
ted vertical line; see Paper-I). This requirement is needed for 
reliable morphological classification and surface photometry. A 
complete and useful log of the observations concerning sky, zero 
point fluctuations, etc. can be found in the CASU website|j. 

All WFCAM data, including proprietary PI data, are 
pipeline-processed by the Cambridge Astronomical Survey 
Unit (CASU) as part of the VISTA Data F l ow System 
(VDFS) development (s ee, Emerson et al.l l2004t llrwin et all 
l2004tlHamblv et al.l2004l for an overview). The data products of 
the CASU pipeline include artifact-corrected images, interleaved 
and stacked as appropriate, confidence maps and catalogues of 
detected objects (both stellar and extended). All of them are 
astrometrically and photometr ically calibrated w ith respect to 
2MASS. We refer the reader to llrwin et al] (120081) for a detailed 
description of the pipeline and iHodgkin et al.N2009l) for a de- 
tailed description and discussion of the photometric calibration. 
An overview of the pipeline and photometric calibration can also 
be found in the WFCAM section of the CASU web siteQ, which 
contains other quality control plots too, showing, for each night, 
the variations of the photometric zero -points, the image see ing 
(FWHM) and sky background level (see iRiello & Irwinll2008l for 
an overview) that proved very useful for assessing the quality of 
the data presented in this paper. 



1 Please refer to WINGS Website for updated details on the survey 
and its products, http : //web . oapd. inaf. it/wings 

2 Note that WFCAM is equipped with a broad K-band filter, not with 
the short-band K s . 



http : //casu . ast . cam . ac . uk/surveys-pro j ects/wf cam/data- 
processing 

4 The web page http://www.ast.cam.ac.uk/~wfcam provides 
access to the WFCAM raw data and includes detailed description and 
useful summaries of the observing and reduction processes. 
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Table 1. Summary of the observations, taken with WFCAM on UKIRT, discussed in this paper. 



Semester 


PATT Ref . 


Alloc. 


Completion 


Observation Dates 


Comments 


2005A 


U/05A/2 


40h 


67% 


Apr29-May 4 
May 14-15(Pq) 
May 28-29(Pq), 2005 


Compromised by poor 

weather, seeing and 
camera misalignment. 


2005 B 


U/05B/12 


40h 


18% 


Sept 30 
Nov 27, 2005 
Jan 6, 2006 (all Pq) 


Severely 
impacted by 
poor weather. 


2006A 


U/06A/12 


lOh 


100% 


May 20&26(Pq) 

June 2(Pq) 
June21&23, 2006 


Good 
observing 
conditions. 










Dec. 14, 2006 




2006B 


U/06B/12 


14h 


65% 


March 20, 28, 31 
April 3-5, 10-11, 2007 (all Pq) 


Average conditions. 


Note "(Pq)" refers to observations taken in the PATT queue, and not by a WINGS project observer. 
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Fig. 2. The central region (2.5'x 2.2', north is up, east is left) of the cluster A2124 K-band mosaic. 



2.1. Pipeline data products 

WFCAM observations are organised in groups, where a group is 
composed of all the images taken with a given filter, dither and, 
optionally, micro-step sequence. For all the clusters presented 
in this paper, observations were taken using a 9-point dithering 
pattern with a 2x2 micro-step sequenc e to partially recover the 
spatial resolution^ dCasali et alj|2007l) . Micro-stepping is done 
by shifting the telescope by n + \ pixels (for the 2x2 pattern) 
and interleaving consists of an algorithm that creates an output 
image that is a regular interwoven pattern of all the input pix- 



3 Micro-stepping and interlacing allow scales of 1/2 and 1/3 of the 
original pixel size, avoiding the PSF under-sampling when good seeing 
conditions occur. In all our observations the effective pixel size is 072. 



els. Some relevant caveats about interleaving for this paper are: 
i) interleaving does not eliminate bad pixels, ii) the PSF often 
varies on short enough timescales to lead to unusually "spiky" 
interleaved PSFs that required additional treatment to deal with. 
The first issue was addressed by adopting a dithering strategy 
that allows an efficient rejection of cosmic rays and bad pixels 
at the stacking phase. We will discuss the second issue in the 
following section. 

Each group is therefore composed of 36 independent Multi- 
Extension FITS files (MEF) with four extensions (one per de- 
tector). Given that the WFCAM four detectors are separated, in 
both directions, by a 90% of their size, it is necessary to take 
four images {i.e., four MEFs) with appropriate offset in order to 
survey a contiguous area of about 0.79 deg 2 (a tile, see Fig|4]i. 
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Fig. 3. The distribution of the median FWHMs (angular vs. physi- 
cal) of all the J- (open blue circle; blue dashed line) and K-band 
(filled red circles; red full line) mosaics: black dotted lines are 
the total J and K distributions. The error bars show the RMS 
of the median seeing variation of each section of the mosaic, as 
shown in FigH] Note that a large error bar means that, due to 
non-photometric nights and/or technical problems with the cam- 
era, the FWHM was unstable among the different exposures mak- 
ing up the final mosaic. The dotted vertical line is the limit for 
quality requirements of WINGS survey spatial physical resolu- 
tion, needed for reliable morphological classification and surface 
photometry. 

Depending on the cluster and the filter, one or more tiles were 
observed to achieve the desired magnitude limit. 

For each image group the pipeline produces a stacked image 
from the interleaved (C'2 sampled) frames coming from each 
dithering point. Each stacked frame is astrometrically and photo- 
metrically calibrated and has an accompanying object catalogue 
and confidence map generated for it. These data products repre- 
sent the input for the final stage of the processing which involves: 

1 . Stacking, for each of the four tile positions, all the individ- 
ual pipeline stacks (typically 2 for J-band observations and 4 
for K-band ones). To avoid confusion we will refer to these 
images as the (four) final stacks. 

2. Mosaicking the four final stacks into a single image covering 
-0.79 deg 2 . 

2.2. Final Stacks and Mosaics 

The final stacks are produced using a stand-alone version of the 
image stacking engine used by the CASU pipeline. The main 
difference of the stand-alone version is that it uses the object 
catalogues for further refinement of the WCS (World Coordinate 
System) offsets that are stored in the FITS headers. In particular: 

1 . The alignment of the input images is driven by the WCS of 
the input frames and then is further refined using the associ- 
ated object catalogues. The first FITS image in the stacking 
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Fig. 4. Due to the WFCAM configuration and the way the obser- 
vations are carried out, the final mosaics have global PSF vari- 
ations in a chessboard pattern fashion. The 'tile' described in 
the text, is composed of 4 exposures (i.e., MEFs), of 4 detec- 
tors each, which gives the 16 sections of the mosaic shown here 
(pixel size is (X'2). The best and the worst global PSF of our mo- 
saics is presented in the top and the bottom panel, respectively. 
For each of the 16 sections, the median FWHM (in pixels) of well 
classified and non saturated stars is reported. 



list is used as a reference and the other images are resampled 
(using nearest-neighbour interpolation) onto the WCS of the 
first image. 

2. All the images are then scaled at the detector level, using ad- 
ditive corrections, to match the background of the common 
overlap area. 

3. The next step involves rejecting bad pixels using the input 
confidence maps and an iterative A:-sigma clipping that uses 
a sky background-calibrated Poisson noise model for the 
frame (which uses a robust MAD0-based sky noise estimator 
to define the equivalent RMS noise level). 

4. Finally, pixels are combined using inverse variance weight- 
ing derived from the combination of the input confidence 
maps and the average noise properties of the image with 
an additional weighting according to the seeing (i.e. weight 
cx seeing^ 2 ). 

A two-dimensional non-linear iterative filter was then ap- 
plied to remove the background from the final stacks. This addi- 
tional step is required because the stacked frames produced by 
the pipeline, although corrected for residual reset-anomaly, gain 
and low-level sky variations, still show background variations 
(o~bg = 0.7 -r 1 .2 A DUs). Finally we used a variant of the drizzle 
algorithm (see e. s. . iFruchter & Ho ok 2002) to remap the origi- 



MAD: median absolute deviation from the median. 
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nal 074 pixels that were interleaved into a finer 072 grid to obtain 
smoother PSFs (see Sec. 12. It . 

The final stacks are eventually combined together to form a 
mosaic, again using a stand-alone version of the pipeline soft- 
ware. The software uses the photometric zero points associated 
with each of the four input images to match them photometri- 
cally and then project them on the same grid using a drizzle-like 
algorithm. To reduce the artifacts caused by the poorer image 
edges, all the input pixels with an associated confidence lower 
than 50% were rejected from the final mosaic. Finally the mosaic 
zero-points are re-measured with respect to 2M ASS using the 
same software used by the CASU pipeline (see iHodgkin et al.l 
120091) . and a final confidence map for the mosaic is produced 
from the confidence maps of the single stacks. 

3. The Photometric Catalogs 

After preparing the final photometrically and astrometrically cal- 
ibrated co-added mosaics (hereafter simply called mosaics) we 
used a custom pipeline of miscellaneous scripts to generate the 
source lists. We will us e SExtractor's MAG-AU TO (with the de- 
fault input parameters, iBertin & Arnoutslll996i see manual for 
details) as the default magnitude throughout this paper, unless 
explicitly otherwise mentioned. 

Hereafter we present a step-by-step schematic description of 
the catalogs production: 

- selection of a polygon enclosing the reliable region of the 
mosaic, whose total area is reported in Tabj3] In practice, the 
edges of the mosaic are excluded. We also exclude objects 
that are dissected by, or are not far enough away from the 
edges to ensure their photometry is unaffected; 

- running the Terapix software WeightWatchers to generate a 
flag mask, to locate the detections inside the good region 
and those with confidence < 70% as deduced from the con- 
fidence maps of the final mosaics: we remind the reader that 
pixels from the single stacks with confidence levels below 
50% were not included in the final mosaic (see previous sec- 
tion); 

- preliminary running of SExtractor, to evaluate the FWHM 
of non saturated stars (CLASS_STAR> 0.95), to choose the 
detection thresholds, and to quickly check the quality of the 
mosaics (astrometry, photometry, number counts ...), as ex- 
plained in section [XT] 

- extraction of several background stamps from the mosaic, 
useful to prepare a background image for photometric errors 
estimates with simulations. A collection of synthetic stars 
and galaxies (30% of exponential disks) is separately added 
to the background image, in an attempt to best reproduce the 
FWHM distribution of the real image. Detection and classifica- 
tion rates of stars and galaxies are computed separately (see 
Tab. [3} and are also used to fine tune the SExtractor input 
parameters (see section [XTT i; 

- final running of SExtractor with the adjusted input param- 
eters and partitioning of the main output catalog into the 
catalogs of stars, galaxies and unknown objects. The lat- 
ter step is achieved here just relying upon SExtractor's 
CLASS_STAR, which ranges from (galaxies) to 1 (stars), in 
the following wajQ: 

7 In paper-II the galaxies value was 0.2: the difference is due to the 
intrinsic differences between near-infrared and optical images. These 
values are estimated using simulations which reproduce the peculiar 
characteristics of the mosaics. 



Stars CLASS_STAR > 0.8 

Galaxies CLASS_STAR < 0.35 
Unknown 0.35 < CLASS.STAR < 0.8; 

- since the CLASS_STAR parameter is not reliable enough in all 
circumstances, in particular when the mosaic is affected by 
strong seeing variations, a further step is necessary in order 
to improve the assignment of each object to the right cata- 
log and to remove from all of them the spurious detections. 
For this task we used a custom interactive tool which pro- 
duces various diagnostic plots, and allows us to distinguish 
between stars and galaxies; 

- final "visual" (interactive) cleaning of the mosaic to correct 
any remaining blatant error both in detection and classifi- 
cation. In this phase the mosaic is displayed with the corre- 
sponding markers of stars and galaxies and, by visual inspec- 
tion, the corrections are made directly on the image, saved on 
disk and then applied to the catalogs. 

In the following sections we will describe in more detail 
some of the aforementioned steps taken as part of the catalogs 
production. 

3.1. Source Detection and Star/Galaxy Classification 

The source detection was performed by running SExtractor 
on the final mosaics, convolved by 2D gaussian filters of size 
chosen to be equal to the median FWHM. In most cases, the low 
level of the background rms obtained with the drizzling tech- 
niques of the CASU pipeline, allowed us to use a 1.5<x clipping 
and a minimum area of 20 adjacent pixels as threshold param- 
eters. These parameter values allowed us to simultaneously ob- 
tain a small number of spurious detections of stars and galaxies, 
and deep enough magnitude limits. In general, for each image, 
the right combination of threshold, minimum area and filter size 
was chosen relying upon the expected number counts of galax- 
ies/stars and on serendipitous visual inspection of marked de- 
tections with the SAO-DS9 display tool. Spurious detections are 
typically misclassified as galaxies by SExtractor, and usually 
result from local background fluctuations and spikes or cross- 
talk of saturated stars. 

The preliminary star/galaxy classification was done relying 
upon SExtractor's stellarity index (CLASS_STAR). The median 
FWHM of the image to be processed (SEEING_FWHM) is the chief 
parameter affecting the SExtractor's star/galaxy separation al- 
gorithm. The sigma-clipping and filtering detection parameters 
used in SExtractor are also quite important, since background 
fluctuations and PSF distortions near the detection threshold can 
introduce uncertainties correlated with such quantities. 

Aiming to test on our mosaics both the detection capability 
of SExtractor and the reliability of its star/galaxy classifier, 
we have produced for each cluster a mock image of simulated 
stars and galaxies (a mix of 70% spheroidal r 1/4 and of 30% 
exponential disks) with a sample background coming from the 
real image, and trying to match as far as possible the distribu- 
tion of the stellar FWHM. With ARTDATA package, synthetic stars 
were modeled using Moffat profiles with /? = 2.5. For spheroidal 
De Vaucouleurs and disk profiles, intrinsic half-light radii in the 
range 1-5 kpc and an Euclidean power law (n=0.3) luminosity 
fun ction were adopted. 

iFoucaud et al.l d2007l) claim that a 20% increase of the value 
of SEEING _JWHM parameter had to be used to correctly match the 
synthetic stars added to the images, and this indicates the diffi- 
culty SExtractor has to correctly classify stars, at fainter mag- 
nitudes, even in good seeing conditions. We further investigated 
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r=FWHM 



AFWHN 



Fig. 5. Distribution of the quantity r = FWHM M in-AFWHM of Eq.© 



for all the near-infrared mosaics. FWHM, 



max/min 



are the maximum 



and the minimum values of the median FWHM, calculated on the 
single detectors of the final stacked tiles, as shown in Figj4] For 
all clusters with r < 0.25 (the region delimited by the vertical 
dotted line) both the star/galaxy classification and the cleaning 
of the catalogs are quite straightforward. 



this effect and found that, while the FWHMs of stars produced by 
the ARTDATA-IRAF package are consistent with the input values, 
their SExtractor estimates ar e ^15% higher . This drawback 
can be dealt with following the iFoucaud et al.l d2007l) prescrip- 
tions of generating mock images with correspondingly dimin- 
ished input values of the FWHM. The results of these simulations 
are presented in Sec. 14.11 where we discuss the completeness of 
the catalogs. 

The PSF spatial variation exemplified in Fig]4] is by far 
the most important effect which might complicate the correct 
choice of CLASS.STAR. In fact, while the resulting CLASS.STAR 
is highly sensitive to the SEEING_FWHM keyword, up to the last 
release of SExtractor has this keyword fixed for the whole 
frame. This means that, in case of strong space variations of the 
FWHM, SExtractor will tend to overestimate the number of stars 
(galaxies) in the mosaic regions where the local FWHM is greater 
(lower) than the median value used for source extraction. 

In Fig|5]we present, for all our clusters in both bands, the 
distribution of the quantity: 



FWHM« 



(FWHM N 



FWHMmin) 



(1) 



where FWHM_max/min are the maximum and the minimum me- 
dian FWHM of the single stacked detectors, respectively, as deter- 
mined from grids like those in Fig fj] The value of r is a qual- 
ity and stability indicator of the PSF for complete mosaics and 
all the observations collected during a given night. Fig|5] shows 
that the majority of our clusters have relatively low values (high 
quality) of this quantity (r <0.25; marked by the vertical dotted 
line). In clusters with r > 0.25 the choice of the SEEING_FWHM 
keyword and other parameters is very difficult, due to many mis- 
classifications being found, and the use of the interactive clean- 
ing procedure becomes fundamental to generate reliable stars 
and galaxies catalogs (see Fig(6]l. 



A1069 



A1795 



/ ■■■•y:.';-v.>v.o-c>./ 



14 16 

J (mag) 




Fig. 6. Plot of the difference between two aperture magnitudes 
(1 arcsec - 3 arcsec diameter) vs. total magnitude for galax- 
ies (blue bold dots) and stars (red light dots) as classified by 
SExtractor before any interactive cleaning. In this diagram the 
star and galaxy loci turn out to be well apart down to J~19.0. 
In the good cases (top panel: A1069, with quite stable seeing 
during observations, t = 0.21), only a few galaxies are misclas- 
sified, and even at bright magnitudes. Alternatively, in cases of 
strong seeing variations (bottom panel: A1795, r = 0.32), many 
more misclassifications of galaxies may be found (bold dots in 
the locus of stars). 



3.2. Interactive Cleaning 

Spurious detections are most frequently found along the over- 
lapping regions of the detectors, along the mosaic edges (which 
are excluded by the initial polygon mask) and in (usually lim- 
ited) regions of high background fluctuations. Cross-talk due to 
saturated stars is another cause of spurious objects. These false 
detections can be found at fixed distances (in symmetric posi- 
tions) from the saturated objects along the read-out direction of 
the detector, and hav e a characteristic "doughnut" shape (see, 
i.e., iDve et al1l2006l) . It was not possible to safely detect and 
delete them automatically, as their number and occurrence is not 
simply correlated with the peak intensity and the area of the sat- 
urated object, and in the final mosaic they can be found either 
in the X or in the Y direction. On the other hand, by exploit- 
ing the SExtractor classification capabilities, it was found that 
only the brightest ones are misclassified as galaxies, and most of 
them are easily deleted by the interactive cleaning procedure de- 
scribed below. Since the automated pipeline produces a separate 
catalog of the saturated stars, during the final visual check, we 
highlight them, allowing easy identification and deletion from 
the final catalogs of their associated cross-talk false detections. 
A negligible fraction of these objects, depending on the number 
of saturated stars in the mosaic, may remain after this last check, 
but are mainly classified as "unknown" sources. 

Since the number of objects is too large for any individ- 
ual analysis, we purged spurious detections and star/galaxy mis- 
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classifications with our interactive tool, which generates several 
plots of different combinations of parameters, like: 

- AWxVs. J,K , where ju max is the surface brightness of the 
brightest pixel of an object and J,K is the total apparent mag- 
nitude MAG_AUTO. 

- J,K vs. log 10 (Isophotal Area), where Isophotal Area is the 
area in pixels of each object at the threshold level. 

- log(FWHM) vs. log l0 (Isophotal Area). 

- J,K( d < l")-J,K(d < 3") vs. J,K, where d is the aperture 
diameter. 

In all these diagrams, the stars populate a rather narrow and 
well defined region, while galaxies are more spread throughout 
the plane (see also, Paper-II). As an example, in Figj6] a plot 
of the difference between two aperture magnitudes (1 arcsec - 
3 arcsec diameter) vs. total J magnitude for galaxies (blue bold 
dots) and stars (red light dots) before any interactive cleaning, is 
presented. It is apparent that for the cluster A1795, due to seeing 
variations and background conditions, a considerable amount of 
objects with CLASS_STAR typical of galaxies actually populate 
the region of stars. In the case of A1069, the more stable PSF 
causes a better classification. Still, at faint magnitudes, a non- 
negligible amount of misclassifications is left. Thus, a further 
step turns out to be necessary in any case to clean the catalogs. 
To this aim, the outliers and/or misclassifications highlighted 
by these diagrams are interactively selected, and, if wanted, the 
cleaning pipeline shows a tile-mosaic of these objects for vi- 
sual inspection to easily select spurious detections. Then, sim- 
ple commands allow us to look at their location in the original 
mosaics, list their parameters from the catalog, delete them or 
change their classification. 

At the end of this process the degree of misclassification of 
relatively bright sources (J< 18.5 and K< 18.0) is less than 1%, 
that is practically negligible (see Sec. 14. U . Of course, going 
to fainter and smaller objects the regions occupied by stars and 
galaxies start to mix up, making the classification more and more 
unreliable. 

The published catalogs will be regularly updated to correct 
for any newly found spurious objects and/or misclassifications, 
therefore users are encouraged to rely upon the latest available 
version of the catalogs. 

3.3. Catalogs Description 

In Tab|2] we present an example of the entries in the near- 
infrared photometric catalogs. The parameters stored for each 
object are the following (in parenthesis we give the name of the 
SExtractor's output parameter): 

- ID: objects internal identification, it is unique for all catalogs 
of the WINGS survey. 

- (aBary^Ban)' equatorial coordinates (J2000.0) of the 
bary center. 

- ((Xpeak,8peak)'- equatorial coordinates (J2000.0) of the bright- 
est pixel. 

- Area: area at the detection threshold. 

- rKnm ■ Kron radius used to compute the MAG_AUT0 magnitude. 

- FWHM: full width at half maximum assuming a gaussian core 
as calculated by SExtractor. 

- b/a : axis ratio of the source. 

- PA : position angle of the major axis (North=0°, measured 
counter-clockwise). 

- Umax '■ Surface brightness of the brightest pixel. 

- MAGJ3EST : SExtractor's best total magnitude estimate. 



- MAGJUJTO : SExtractor's Kron (total) aperture magnitude. 

- MAG(4.3kpc), MAG(10.8kpc) , MAG (21 . 5kpc): magni- 
tudes within apertures of diameter 4.3kpc, 10.8kpc and 
21.5kpc, respectively, measured at the clusters' redshift. 
Since we will adopt in future papers of the survey Hq = 70 
(at variance with Paper-II, where Hq = 75 has been used), 
the apertures in kpc are slightly different from those given 
in Paper-II (4/10/20 kpc): however, the value in arcseconds 
remains the same. 

- HkG(&.'7), MAG(l'.'Q), MAG(1'.'6) . . . : magnitudes mea- 
sured within 8 different fixed aperture^] in arcsec. 

- FLAG_SEX : SExtractor's FLAG keyword. 

- CLASS_STAR : SExtractor's stellarity index. 

- WINGS_FLAG: summarizing flag column reporting pipeline 
classification and photometric quality of the objects, using 
the following prescription: 

WINGS _FLAG = ai + 2a 2 + 4a 3 + 8a 4 + 16a 5 + 32a 6 
a i — 1 if classified as star 
«2 = 1 if classified as galaxy 
«3 = 1 if classified as unknown 
«4 = 1 if weakly affected by neighbouring halo 
«5 = 1 if strongly affected by neighbouring halo 
«6 = 1 if in an area of confidence <70% 
where the classification is done relying upon our custom in- 
teractive tools, the neighbouring halo can be due to satu- 
rated stars or nearby extended objects, and the confidence 
level is extracted from the final confidence maps generated 
by CASU pipeline (see section [2~2| >. 

The complete/updated list of stored parameters, together 
with SExtractor configuration files, are available on the 
WINGS website http://web.oapd.inaf.it/wings/, while 
the exact definition of each parameter can be found in the 
SExtractor manual dBertin & Arnoutsl 19961) . The photometric 
catalogs, in the form of Tabf2]are available in electronic format. 
Note that the same caveats and/or qualifications given in Paper-II 
hold for these photometric and geometric parameters. 

The final version of all WINGS survey catalogs will be a 
comprehensive cross-matched source-list, allowing multiple cri- 
teria queries with Web based tools and Euro-VO facilities, and 
simple identification of objects which have measurements and 
calculated quantities belonging to different branches of the sur- 
vey (U,B,V,J,K photometry, spectroscopy). 

4. Data Quality 

As already mentioned in previous sections, the bulk of the re- 
duction has been performed at CASU with the custom designed 
pipelines for UKIRT-WFCAM. Hereafter we present the results 
of the quality tests performed by our pipeline on source extrac- 
tion and classification, both in real and simulated mosaics. 

4.1. Completeness 

FigfT] illustrates the global detection and successful classifi- 
cation rates (for both stars and galaxies) evaluated running 
SExtractor on the mock images, as explained in the previous 
Sec. l3.1l (see also Tabf3]l. Full and dashed lines are the fraction of 

8 Three of these fixed apertures have been chosen to match the multi- 
fiber spectroscopy of the WINGS fields. 1'.'6 is the projected diameter of 
the fibers using the Autofib2@WHT, while using the 2dF@ A AT results 
in a diameter that varies radially in the field from 2'.' 16 in the center to 
2'.'0 at the edges. 
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Table 2. Example of entries in the on-line version of J and K photometric catalogs. 



ID 




<5flary 


"Peak 


&Peak 


Area 


r Kron 


FWHM 


b/a 


PA 


rmax 


(deg) 


(deg) 


(deg) 


(deg) 


(arcsec 2 ) 


(") 


(") 




(deg) 


(mag arcsec -2 ) 


WINGSJ 1 03934.24-082406.9 


159.89270 


-8.40195 


159.89267 


-8.40193 


8.08 


0.70 


1.90 


0.73 


99.34 


18.70 


WINGS J 1 03944. 1 9-0824 1 0.8 


159.93416 


-8.40298 


159.93414 


-8.40300 


4.68 


0.70 


1.65 


0.87 


10.52 


19.36 



MAGJ3EST 


MAG_AUT0 


MAG(4.3kpc) 


MAG(10.8kpO 


MAG(21.5kpc) 


MAG(S'.'7) 


BAGCl'.'S) 


BAG(1'.'6) 


HAG(2'.'8) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag ) 


(mag) 


17.63 


17.63 


17.69 


17.62 


17.84 


19.87 


19.17 


18.41 


18.13 


18.34 


18.34 


18.36 


18.35 


18.03 


20.50 


19.82 


19.07 


18.80 



MAG(2'.'16) 

(mag) 


HAG(3'.'8) 

(mag) 


MAG(4'.'S) 
(mag) 


MAG(576) 

(mag) 


FLAG.SEX 


CLASS.STAR 


WINGSJLAG 


18.05 


17.75 


17.64 


17.60 





0.02 


2 


18.71 


18.44 


18.31 


18.22 





0.01 


2 



The full table, for each cluster and filter, is available in electronic form at the CDS. 
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Fig. 7. The global average detection (full lines) and success- 
ful classification (dashed lines) rates, obtained from simulations 
for the WINGS-NIR survey in the J- (left panels) and K-band 
(right panels). The dotted lines represent the 95, 90 and 75 
percent of completeness. The dotted-dashed line in the upper 
panels represent the fraction of stars misclassified as galaxies 
by SExtractor. The bottom panels show that, even at bright 
magnitudes, some galaxies (a small fraction, indeed) are miss- 
classified as stars. This is likely due to the wide range of effective 
radii adopted in the simulations, which can sometime produce 
unrealistically compact luminosity profiles. In real mosaics, such 
misclassifications are easily corrected with our interactive tool. 



the total input objects, detected and correctly classified, respec- 
tively. The 90% level of successful classification is reached at 
J=19.0 and K=18.2 for stars, while the limits become J=19. 5 and 
K=18.5 for extended sources. Moreover, the 90% detection rate 
limits for galaxies is reached at J=20.5 and K=19.4. In Tabj3]we 
report the detection and classification 90% completeness limits 
for both stars and galaxies for each cluster and band, as extracted 
from our simulations. 



It is clear that SExtractor is able to classify correctly al- 
most all the galaxies detected in simulations (apart from rare 
and very compact objects), while this is not true for point-like 
sources. This is likely due to the seeing variations across the mo- 
saics (which in a few cases is relevant, see Figs |4|5t . In FigJTJ the 
three dotted horizontal lines mark the 95%, 90% and 75% com- 
pleteness limits of our survey, while the red dashed-dotted line 
in the upper panels represent the fraction of stars misclassified as 
galaxies when relying only upon SExtractor's CLASS.STAR. It 
can be seen that for J,K>18.0 the frequency of such misclassifi- 
cation becomes relevant (see also Fig(6]) and rapidly increases. 
However, we have to consider that, for faint magnitudes, the 
galaxy population becomes gradually dominant (field galaxies) 
and the contribution to the galaxy counts coming from the (usu- 
ally small) fraction of misclassified stars should be in any case 
negligible (see, for e.g. jBerta et alj[2006b . 

An empirical (a posteriori) check of completeness can be 
performed, for each cluster in each band, comparing the num- 
ber counts of classified objects in our catalogs, with the expected 
number counts of field galaxies and stars. In Fig[8]we show rep- 
resentative cases for J- (left panels) and K-band (right panels). 
The red dashed lines are the expected number counts of stars 
in the area of the c luster, calculated with the TRILEGAL code 
dGirardi et al.ll2005l) in the photometric system of WFCAM. The 
black solid lines are the number counts of field galaxies taken 
from the UKIDSS Ul tra Deep Survey (see. iHartlev et al.ll2008h 
lLawrence et al. 1 12007b . The open red and the full blue dots are 
the WINGS number counts for stars and galaxies, respectively, 
as calculated from our source lists, with their Poissonian error 
bars. In general, the agreement between the expected star counts 
and our catalogs is outstanding, confirming the excellent per- 
formance of the TRILEGAL code as a model of the Galaxy. 
However, for 6 clusters out of 28, we noticed an excess in the 
real star counts at magnitudes brighter than 18 (see for instance 
A1644 in Fig[8]). Since all these clusters are found to lay approx- 
imately in the direction of the center of the Galaxy (albeit with 
\b\ > 20"), where a higher fraction of halo stars is encountered, 
we assume that the model might slightly break down in this di- 
rection. Our galaxy counts agree with the field number counts 
only at faint magnitudes, as expected, showing the presence of 
cluster galaxies at brighter magnitudes. 

The excess of number counts shown in Figj8] revealing the 
cluster members, can also be seen in Fig|9] where the number 
counts of all the detected sources (full line), is compared with the 
cumulative (stars+field galaxies) counts expected for the area of 
the clusters (dashed line). The approximate magnitude where the 
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Fig. 8. Number counts in J- (left panels) and K-band (right pan- 
els) compared with TRILEGAL model s of the Galactic stel- 
lar distribution from lGirardi et al.l d2005l) ( red dashed lin e s) and 
field galaxy counts from UKIDSS-UDS lHartley et al.l (l2008l) 
(solid black lines). The blue full and the red open dots are the 
galaxy and stars number counts of the corresponding WINGS 
catalogs. Poissonian error bars are mostly of the dimension of 
the symbols. 



turnover in the counts occur is a good estimate of the observed 
completeness limit, and is reported in Tab|3] together with the 
theoretical magnitude limit (dotted-dashed vertical lines) calcu- 
lated with the following formula: 

mi im = ZPT - 2.5 log 10 [v (cr BG + 1) A min ] (2) 

where v is the relative threshold cut in units of background rms 
(o~bg), A m ; n is the minimum number of contiguous pixels re- 
quired for detection and ZPT is the mosaic zero point magnitude, 
normalized to one second exposure time and airmass corrected: 

ZPT = MAGZPT + 2.5 log 1() (5) - (x - 1) • k (3) 

where MAGZPT is the zero-point magnitude keyword found 
in the FITS header of the mosaic, resulting from the calibra- 
tion with 2MASS performed at CASU. The second refers to 
the WFCAM 5s exposure time, the term \ an d k variables are 
the average airmass and extinction (keywords AIRMASS and 
EXTINCT in FITS file), respectively. 

This theoretical limit corresponds to the magnitude of an ob- 
ject consisting of A m ; n contiguous pixels with ADU counts of 
v{ctbg + 1); an d gives an idea of the overall depth of the mo- 
saic imaging, since it links photometric and detection properties 
together. The dash-dotted vertical lines in Fig|9]show that these 
limiting magnitude values are quite consistent with the faintest 
detected objects. 

At the final step, the consistency among the previous three 
diagrams (Figs. |9]|7] and [8]l is used as an a posteriori quality 
check, as well as as a good empirical tool for possible refine- 
ment of SExtractor's input parameters. Usually, when this 



Fig. 9. Magnitude counts of all detections in J- (left panels) and 
K-band (right panels), for the deepest (upper panels) and the 
shallowest (bottom panels) cluster mosaics. The dash-dotted ver- 
tical lines correspond to the estimated theoretical magnitude lim- 
its from Eq.©; the turnovers in the counts given in this figure in- 
dicate the detection completeness for these cluster images. Both 
values are reported in TabJUfor all clusters. The dashed line is 
the cumulative stars+field galaxies distribution expected to be 
found in the area of the mosaic. 



ideal combination is achieved, objects are detected down to the 
completeness limit without populating the source lists with un- 
wanted spurious detections. 

4.2. Astrometry 

Most of the non-linear distortion of WFCAM over the entire 
field is accounted for, by a cubic radial term in the astromet- 
ric solution. The CASU pipeline processes the raw images con- 
sidering the differential field distortion too, giving at the end an 
astrometric error usually below 50maj^ ( for further details re- 
fer to CASU website, and lDve et afll200l5l) . We are able to reach 
that precision only in relative astrometry, i.e. difference of co- 
ordinates in two bands of the same field. After the coadding 
of all the MEFs into a single mosaic, an additional astromet- 
ric and photometric re-calibration check with point-like sources 
from the 2MASS catalogue is performed. Fig[lO]is a visualiza- 
tion of the astrometric precision and accuracy for our entire clus- 
ter collection: upper panels show the overall spread of the right 
ascension and declination differences between the positions of 
common point-like sources in our catalogs and in 2MASS (left 
panel) or UCAC2 (right panel). Only sources with a photomet- 
ric error lower than 0. 1 mag were considered. The overall zero 
point shift is negligible for all our applications, and the rms is 
consistent with UKIRT- WFCAM standard requirements, being 
of the order of lOOmas (RMS=1 12mas). 
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Fig. 10. Astrometry of the WINGS-NIR survey, compared with 
2MASS point-source and UCAC2 catalogs (upper panels) and 
histograms of the distances (bottom panels). The overall preci- 
sion is of the order of lOOmas and nearly 70% of the offsets 
calculated are below one pixel distance. The high precision as- 
trometry of our survey allows safe cross-matching of catalogues 
extracted from other surveys. 

While the original stacks of each tile use a Zenithal 
Polynomial projection (ZPN) astrometric solution, the final 
coadded mosaics are expressed in the standard gnomonic tan- 
gential projection (TAN). Thanks to the accurate and precise as- 
trometric solutions applied to our mosaics, cross-matched source 
lists in different bands (optical and U bands included) are very 
easily obtained, and the fraction of mismatches is negligible. 

4.3. Photometry 

In Tab|4]a summary of the properties of the catalogs is given, and 
the surface brightness limits, calculated in the following way, are 
reported: 

Wim = ZPT - 2.5 log 10 [v <t B g] + 2.5 log 10 [A"] (4) 

where A" is the pixel area in arcseconds. This relation gives the 
minimum surface brightness a pixel can have due to the sigma 
clipping chosen for the specific mosaic. The output of Eq.© 
is obviously changing from cluster to cluster, but it generally 
settles at /j.j\ lm « 22 and //Kiim ~ 21 . 

Photometry at CASU is currently based on 2MASS, via 
colour equations converting 2MASS magnitudes to the WFCAM 
photometric system. The m ost recently released photometric cal- 
ibrations are given by (see. lHodgkin et ai1l2009l) . Neglecting the 
interstellar extinction, they are as follows: 

JwFCAM = J2MASS ~ 0.065(J2MASS ~ H2MASS) (5) 
K-WFCAM = K2MASS + 0.010(J2MASS - K2MASS) (6) 

Due to the improved precision in the fitting algorithm, these 
equations are different from those given in the early-data-release 



Fig. 11. The 2MASS photometry is converted to the 
WINGS/WFCAM photometric system with Eqs.dBToT) and 
then compared with WINGS magnitudes, for the whole col- 
lection of our NIR mosaics, separately in J- (upper panel) and 
K-band (lower panel). For a correct comparison with 2MASS 
we use aperture (3'.'0 diameter) corrected magnitudes in the 
Y axis. The absence of a significant zero point shift and the 
relative tightness of the magnitude sequence demonstrates the 
quality of the photometry in the WINGS survey. 

bv lDve et all (|2()06). Different tests carried out at CASU suggest 
that the 2MASS calibration is delivering photometric zero-points 
at the +2% level, confirming the excellent results achieved by 
the 2MASS survey team in ensuring a reliable all-sky accurate 
calibration. While coadding the single MEFs to obtain the fi- 
nal mosaic, a photometric re-calibration check based on 2MASS 
catalogs is performed again, to assure a spatially homogeneous 
zero point throughout the mosaic. FigQl] shows the difference 
between 2MASS and WINGS aperture corrected magnitudes of 
point-like sources (after applying Eqs|5]and|6]l, vs. WINGS to- 
tal magnitudes in the J (upper panel) and K (lower panel) bands, 
for the complete survey. No significant zero point shift is recov- 
ered, confirming the accuracy of the photometric calibration of 
our WINGS-NIR survey. 

Photometric errors assigned by SExtractor are unrealisti- 
cally small, because they are based only on the photon-noise 
statistics. The use of the so called confidence-maps (weight- 
maps for SExtractor users) helps in computing more real- 
istic values of photometric uncertainties. However, we have 
verified that this latter approach still leads to underestimates 
of the photometric errors. Therefore, we preferred to use the 
simulations to recover the effective precision of our measure- 
ments. Nevertheless, SExtractor errors will be available in the 
WINGS website queries forms. 

In Fig[T2]we present the global results (J and K band to- 
gether) of photometry checks from the simulations. In the top 
panel the differences between input and output SExtractor 
magnitudes as measured from the mock images for stars (red 
starred dots), exponential disks (open magenta circles) and 
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Fig. 12. Global photometric precision and accuracy as obtained 
from simulations for stars (red starred dots), disks (open magenta 
circles) and spheroids (full blue circles). (Top panel) average 
differences between the input and SExtractor magnitudes for 
bins of 0.5mag of input total magnitude. (Bottom panel) rms of 
the previous differences. While the magnitudes of stars and disks 
are perfectly recovered, there is a systematic shift of 0.2mag in 
the spheroidal r 1//4 magnitudes, due to the loss of light in the 
wings of simulated galaxies. 

spheroidal r'^ 4 (blue filled circles) are plotted for 0.5 mag bins 
of total magnitude. The agreement is quite good, apart from 
the constant shift of 0.2mag for spheroidal galaxies even at 
brigh t magnitudes. This effec t is well known from the litera- 
ture dFranceschini et all 1 19981) and is also recalled in Paper-II. 
It is caused by the loss of light from the wings of profiles of 
spheroidal r 1 ^ 4 galaxies in the background noise. The bottom 
panel of the same figure presents the overall precision of our 
photometry (the RMS of measurements is plotted for 0.5mag 
bins, the symbols are as above), with the corresponding poly- 
nomial fit for spheroids (blue full line) and exponential disks 
(magenta dashed line) to be used as the RMS value for all mag- 
nitudes of the catalogs: 



RMS early = -3.675908 + 0.778 lUm ear i y + 
-0.054269»4 rfy + 0.001263m; 
RMSi ate = -3.765820 + 0.806960m tof , + 
-0.057414m^ e + 0.001 



arty 



(7) 
(8) 



where m ear [ y /i ate can be either the total or the aperture magnitude 
of the corresponding object. Due to the way these errors have 
been computed, the statistical Poissonian error can be ignored. 
Moreover Eq|7]is adequate for the errors associated with the "un- 
known" objects category, while Eq|8]can safely be used for point 
like sources. 

FigfT2l shows that the global rms is below 0.1 mag down to 
18.3 mag for stars and disks, and almost 18.0 mag for spheroids. 
It is clear that for J-band mosaics this limit can sometimes reach 
18.8 mag, while for the K-band it is found at «T7.5 mag, de- 
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Fig. 13. Color magnitude diagrams of classified galaxies for a 
selection of WINGS-NIR clusters. The relatively thin red cluster 
sequences confirm the accuracy of the photometry provided by 
SExtractor. The large number of detections obtained at fainter 
magnitudes are the field background galaxies and demonstrate 
the deepness of our survey. There are also indications of red se- 
quences for the field galaxies at bright magnitudes, indicating 
the presence of background clusters at higher redshift. The big 
red full circles are the spe ctroscopically selected members of the 
clusters (ICava et al.ll2Q09h . 



pending on overall quality of the image. Given the overall rates 
of successful classification deduced from Fig|7]and Tabj3] it is 
more than conservative to assign those levels of photometric er- 
rors beyond 18.0 mag and 17.5 mag for the J- and K-band, re- 
spectively. 

As an internal check of photometric accuracy, we show in 
FigE] two representative examples of total MAG_AUT0 color- 
magnitude diagrams of galaxies from our survey. The tightness 
of the red sequences (the RMS of the spectroscopically con- 
firmed members of the clusters is 0.03mag and 0.05mag) demon- 
strates the internal consistency of th e WINGS-NIR photom etry 
(these RMS values are found also bv lEisenhardt et al]|2007l) . 

5. Summary 

In this paper we have presented the first data release of the 
WINGS-NIR survey, comprised of J- and K-band photometric 
catalogs for a sub-sample of 28 nearby galaxy clusters belonging 
to the WINGS optical survey. The detected sources have been 
classified as stars, galaxies and unknown objects, and for each 
of them we give positions, geometrical parameters and differ- 
ent kinds of total and aperture magnitudes in both bands (when 
available). 

Due to the variations of the seeing FWHM over the large ar- 
eas of our mosaics and the poor reliability of SExtractor's 
CLASS_STAR parameter, we used a specifically designed inter- 
active pipeline to improve star/galaxy classification. The near- 
infrared WINGS data consists of nearly one million detected 



12 



T. Valentinuzzi et al.: WINGS III: Near-Infrared Catalogs 



sources, with 150,000 and 500,000 reliably classified stars and 
galaxies, respectively. 

This unique collection of data turns out to be 90% com- 
plete in detection limits at J=20.5 and K=19.4, and in classi- 
fication limits at J=19.0 and K = 18.5 (limits vary from cluster 
to cluster depending on seeing conditions and total integration 
time), values to be compared with 2MASS, J=15.8 and K=14.3 
for point-like sources, UKIDSS Large Area Survey J=20.0 and 
K=18.4, and UKIDSS Ultra Deep Survey J=25.0 and K=23.0 
(for further details, s ee ISkrutskie etail 120061: iDve et all 120061: 
lLawrence etai]|2007l) . 

Using this new, extensive and comprehensive catalog, we 
will study several different properties and characteristics of low- 
redshift cluster galaxies. The near-infrared data presented in this 
paper will be combined with our optical photometry, morphol- 
ogy and spectroscopy catalogs to systematically study the de- 
pendence of these cluster galaxy properties on their stellar mass 
(since this is traced by the J and K bandpasses). Furthermore, 
the study of the J,K luminosity functions (known to be a good 
tracers of the stellar mass function), will allow us to estimate the 
distribution of the stellar mass-to-light ratio as a function of the 
cluster-centric distance, for different morphological types. For 
the galaxies in our sample for which we can derive reliable near- 
infrared surface photometry and structural parameters (e.g., R e , 
< He>), we will investigate the behaviour of our sample with re- 
spect to the various scaling relations (i.e., Fundamental Plane, 
Kormendy relation). Broad-band spectral energy distributions 
will be generated for the galaxies in our sample by combining 
the optical and near-infrared photometry, giving us further infor- 
mation on stellar content and cluster membership. Also, using 
our cluster sample, we will be looking at what effect the use of 
a near-infrared total cluster luminosity has on the Fundamental 
Plane of Galaxy Clusters. 

The final version of the complete WINGS survey cata- 
logs will be a comprehensive cross-matched source-list, allow- 
ing multiple criteria queries with Web based tools and Euro- 
VO facilities. Through the website the identification of objects 
which have measurements and calculated quantities from differ- 
ent branches of the survey (U,B,V,J,K photometry, surface pho- 
tometry, morphology and spectroscopy) will be made easier by 
means of specific web applets. 
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Table 3. WINGS-NIR clusters sample with relevant useful quantities. 



Cluster 


RA 


DEC 


Redshift 


Lx/10 44 


Pixel 


RUN 


FWHM" 


Stacks 


Total Area 




hh:mm:ss 


dd:mm:ss 




erg s -1 


kpc 




ixiin 


max 




. - 


Mpc 2 


A119 


10:39:43.4 


-08:41:12.4 


0.0442 


1.65 


0.174 


05B-K 


0.79 


0.88 


16 


0.781 


7.66 


A376 


02:46:03.9 


+36:54:19.2 


0.0484 


0.71 


0.190 


05B-J 


1.05 


1.17 


8 


0.777 


9.09 


A500 


04:38:52.5 


-22:06:39.0 


0.0670 


0.72 


0.256 


05B-J 
05B-K 


1.28 
1.21 


1.44 
1.42 


8 

16 


0.776 


16.48 


A602 


07 53 26.6 


+29 21 34.5 


0.0619 


0.57 


0.238 


06B-K 


0.88 


0.94 


12 


0.778 


14.29 


A957x 


10:13:38.3 


-00:55:31.3 


0.0436 


0.40 


0.172 


05A-J 
06B-K 


1.02 
0.75 


1.47 
0.92 


8 
8 


0.779 
0.782 


7.47 
7.50 


A970 


10:17:25.7 


-10:41:20.2 


0.0587 


0.77 


0.228 


06B-J 
05A-K 


0.91 
1.06 


1.11 
1.26 


8 

16 


0.777 
0.775 


13.09 
13.05 


A 1069 


10:39:43.4 


-08:41:12.4 


0.0650 


0.48 


0.250 


05A-J 
05A-K 


0.98 
0.88 


1.19 
1.12 


8 
16 


0.783 
0.782 


15.86 


A1291 


11:32:23.2 


+55:58:03.0 


0.0527 


0.22 


0.206 


06B-J 
05A-K 


0.80 
0.99 


0.93 
1.15 


8 

16 


0.778 


10.70 


A1631a 


12:52:52.6 


-15:24:47.8 


0.0462 


0.37 


0.182 


05A-J 
05A-K 


0.89 
0.82 


1.15 
1.10 


8 
16 


0.781 
0.780 


8.38 


A 1644 


12:52:52.6 


-15:24:47.8 


0.0473 


1.80 


0.186 


06B-J 
05A-K 


0.88 
0.89 


0.98 
1.15 


8 

16 


0.778 
0.772 


8.72 
8.65 


A 1795 


13:48:52.5 


+26:35:34.6 


0.0625 


5.67 


0.240 


05A-J 
05A-K 


1.42 
0.95 


1.65 
1.57 


8 
8 


0.780 
0.387 


14.56 
7.22 


A1831 


13:59:15.1 


+27:58:34.5 


0.0615 


0.97 


0.238 


06B-J 
05A-K 


0.72 
0.93 


0.91 
1.07 


12 
16 


0.778 
0.779 


14.23 


A1983 


14:52:55.3 


+ 16:42:10.6 


0.0436 


0.24 


0.172 


06B-J 
05A-K 


0.90 
0.79 


1.03 
0.98 


8 

16 


0.779 
0.772 


7.47 
7.40 


A1991 


14:52:55.3 


+ 16:42:10.6 


0.0587 


0.69 


0.228 


05A-J 


0.88 


1.00 


8 


0.779 


13.12 


A2107 


15:39:39.0 


+21:46:58.0 


0.0412 


0.56 


0.162 


06A-J 
05A-K 


0.95 
0.85 


1.19 
1.03 


8 
16 


0.778 
0.748 


6.62 
6.36 


A2124 


15:44:59.0 


+36:06:33.9 


0.0656 


0.69 


0.252 


05A-J 
06A-K 


1.38 
0.80 


1.18 
1.03 


8 
16 


0.780 
0.782 


16.05 
16.09 


A2149 


16:01.28.1 


+53:56:50.4 


0.0679 


0.42 


0.260 


06A-K 


0.79 


0.98 


20 


0.781 


17.11 


A2169 


16:13:58.1 


+49:11:22.4 


0.0586 


0.23 


0.226 


06A-J 
05A-K 


0.89 
0.83 


1.03 
0.96 


8 
16 


0.780 
0.771 


12.91 
12.76 


A2382 


21:51:55.6 


-15:42:21.3 


0.0618 


0.46 


0.238 


05A-K 


0.89 


1.25 


16 


0.770 


14.13 


A2399 


21:57:01.7 


-07:50:22.0 


0.0579 


0.51 


0.224 


05B-K 


0.91 


1.14 


16 


0.781 


12.70 


A2457 


22:35:40.8 


+01:29:05.9 


0.0594 


0.73 


0.230 


05B-K 


0.84 


1.00 


16 


0.780 


13.37 


A2572a 


23:17:12.0 


+ 18:42:04.7 


0.0403 


0.52 


0.160 


06A-K 


0.79 


0.94 


16 


0.781 


6.48 


A2589 


23:23:57.5 


+ 16:46:38.3 


0.0414 


0.95 


0.164 


06A-K 


0.87 


1.07 


16 


0.781 


6.81 


IIZW108 


21:13:55.9 


+02:33:55.4 


0.0493 


1.12 


0.192 


05B-J 
05A-K 


0.98 
0.87 


1.66 
1.03 


8 

16 


0.777 
0.780 


9.28 
9.32 


MKW3s 


15:21:51.9 


+07:42:32.1 


0.0450 


1.37 


0.176 


HAD T 

06A-K 


0.92 
0.88 


1.03 
1.06 


Q 

o 

16 


0.780 
0.781 


7.83 


RX1022 


10:22:10.0 


+38:31:23.9 


0.0491 


0.18 


0.192 


05B-J 
06B-K 


0.91 
0.80 


1.22 
0.95 


8 
8 


0.779 
0.782 


9.30 
9.34 


RX1740 


17:40:32.1 


+35:38:46.1 


0.0430 


0.26 


0.170 


05A-K 


0.91 


1.23 


12 


0.778 


7.28 


Z8338 


18:11:05.2 


+49:54:33.7 


0.0473 


0.40 


0.186 


06A-J 
05A-K 


0.80 
0.79 


0.95 
0.92 


8 
16 


0.778 


8.72 



RA, DEC: coordinates of the Brightest Cluster Galaxy 
Redshift: from NED. 
Lx: X-ray luminosity 

Pixel: pixel size in kpc at the given redshift. 

RUN: semester and passband of observation. 

FWHM: minimum and maximum estimation as shown in FigO 

Stacks: Interleave Stacks per mosaic used. 

Total Area: the effective area of the mosaic. 
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Table 4. WINGS-NIR catalogs useful parameters. 



Cluster 


DM 


Band 


mii m 


/*tbr6sh 


Cum 


Stars 


Gals 


Notes 






mag 


mag 


mag/" 2 


Det 


Det 


Clas 


Det 


Clas 




A119 


36.46 


K 


21.2 


21.02 


20.0 


19.8 


18.2 


19.3 


19.0 




A376 


36.66 


J 


21.5 


22.11 


20.2 


20.4 


19.0 


20.0 


19.6 


1 


A500 


37.40 


J 
K 


21.6 
20.8 


21.93 
21.02 


21.0 
19.7 


20.0 
19.2 


18.0 
17.7 


19.6 
18.7 


19.2 
18.3 


1,2 
2,4 


A602 


37.22 


K 


21.2 


21.68 


20.0 


20.3 


18.7 


19.3 


19.2 


1,3 


A957x 


36.43 


J 
K 


21.5 
21.1 


22.11 
21.49 


20.3 
20.0 


20.3 
20.0 


18.5 
18.7 


19.8 
19.7 


19.4 
19.4 


2,4 


A970 


37.10 


J 
K 


21.8 
20.9 


22.31 
21.13 


21.0 
19.7 


20.6 
19.5 


18.6 
17.7 


20.0 
18.9 


19.5 
18.5 


2 

1.4 


A 1069 


37.33 


J 
K 


21.9 
20.6 


21.95 
21.32 


21.2 
19.2 


20.2 
20.3 


18.6 
18.7 


19.8 
19.7 


19.5 
19.4 


5,6 


A1291 


36.85 


J 
K 


21.8 
20.6 


22.37 
21.29 


21.0 
19.5 


21.1 
19.5 


19.2 
18.2 


20.3 
19.0 


20.2 
18.8 


1,2,5 


A1631a 


36.56 


J 
K 


21.2 
20.7 


22.09 
21.02 


20.0 
19.1 


20.1 
19.7 


18.5 
18.4 


19.5 
19.3 


19.1 
18.7 


1.3 


A1644 


36.61 


J 
K 


21.8 
20.6 


21.98 
20.90 


20.2 
19.0 


20.7 
19.5 


19.2 
18.1 


20.0 
18.9 


19.6 
18.7 


1,4 
1,4,6 


A1795 


37.24 


J 
K 


22.0 
20.8 


22.27 
21.07 


21.5 
19.0 


20.1 
19.1 


18.0 
18.0 


19.8 
18.7 


18.0 
18.5 


2,5 
2 


A1831 


37.20 


J 
K 


21.7 
20.5 


22.34 
21.19 


21.1 
19.3 


21.0 
19.3 


19.5 
17.5 


20.5 
19.0 


20.1 
18.6 


6 

1,6 


A1983 


36.43 


J 
K 


21.9 
20.2 


22.63 
20.74 


20.6 
18.9 


21.0 
19.0 


19.5 
18.0 


20.3 
18.5 


20.0 
18.5 


1,4 
4 


A1991 


37.10 


J 


21.7 


22.10 


21.1 


20.8 


19.0 


20.0 


19.5 




A2107 


36.30 


J 
K 


21.7 
20.8 


22.22 
21.08 


21.2 
19.2 


20.7 
19.2 


19.1 
18.0 


20.0 
19.3 


19.6 
19.1 


2 

1,4 


A2124 


37.35 


J 
K 


21.7 
21.1 


22.29 
21.36 


21.2 
20.0 


20.2 
19.8 


18.7 
18.3 


20.0 
19.2 


19.7 
19.1 




A2149 


37.43 


K 


20.9 


21.11 


20.0 


19.7 


18.3 


19.3 


19.0 


6 


A2169 


37.09 


J 
K 


21.7 
20.6 


22.14 
20.99 


21.0 
19.7 


20.8 
19.5 


19.1 
18.2 


20.1 
19.3 


19.8 
18.9 




A2382 


37.21 


K 


20.3 


20.57 


19.5 


19.2 


17.7 


18.8 


18.3 


2 


A2399 


37.07 


K 


20.8 


21.04 


19.6 


19.5 


18.2 


19.2 


18.6 


1,4 


A2457 


37.12 


K 


20.9 


21.14 


19.8 


19.7 


18.0 


19.3 


18.7 


4 


A2572a 


36.25 


K 


21.0 


21.25 


19.9 


19.7 


18.4 


19.3 


19.2 


2,5 


A2589 


36.31 


K 


21.1 


21.36 


20.0 


19.8 


18.0 


19.5 


18.8 


2,5 


IIZW108 


36.70 


J 
K 


21.5 
21.03 


21.83 
21.32 


20.4 
19.8 


20.5 
19.8 


19.0 
18.3 


19.8 
19.5 


19.0 
19.0 


1,2,4 
1,4,6 


MKW3s 


36.50 


J 
K 


21.9 
21.0 


22.56 
21.35 


20.5 
20.0 


21.0 
19.7 


19.6 
18.3 


20.3 
19.3 


20.0 
18.7 


1 

2,4 


RX1022 


36.69 


J 
K 


21.5 
21.2 


22.10 
21.53 


20.2 
20.1 


20.5 
19.8 


18.9 
18.2 


19.9 
19.3 


19.6 
19.2 


2 
2 


RX1740 


36.40 


K 


20.6 


20.86 


19.1 


19.1 


17.8 


18.7 


18.2 


1,2,4 


Z8338 


36.61 


J 
K 


21.7 
20.7 


22.16 
21.06 


20.9 
19.5 


20.7 
19.6 


19.2 
18.6 


20.0 
19.3 


19.9 
18.9 


2 



The Cluster column identifies the cluster, the DM column is the cluster distance modulus in magnitudes, the mi im column is the theoretical 
magnitude detection limit calculated with Eq.10, yUthrcsh column is the minimum surface brightness at the detection threshold cut calculate with 
Eq.©, the cumulative detection column ("Cum Det") lists the turnover magnitude of the number counts diagram for all the detected sources (see 
FigHJ, "Stars" and "Gals" columns represent the average 90% completeness detection and classification limits for stars and galaxies as deduced 
from simulations. 
NOTES LEGEND: 

1:= patchy and/or noisy background, 2:= strongly variable PSF, 3:= many spurious detections 
4: = excess of classified stars, 5:= deficient in classified stars, 6:= exceptionally bright star(s) 
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